Line-Defect Patterns of Unstable Spiral Waves in Cardiac Tissue 
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Spiral wave propagation in period-2 excitable media is often accompanied by line-defects, the locus of points 
with period- 1 oscillations. Here we investigate spiral line-defects in cardiac tissue where period-2 behavior has a 
known arrhythmogenic role. We find that the number of line defects, which is constrained to be an odd integer, 
is three for a freely rotating spiral, with and without meander, but one for a spiral anchored around a fixed 
heterogeneity. We interpret analytically this finding using a simple theory where spiral wave unstable modes 
with different numbers of line-defects correspond to quantized solutions of a Helmholtz equation. Furthermore, 
the slow inward rotation of spiral line-defects is described in different regimes. 



Spiral waves are observed in extremely diverse physical and 
biological excitable media and are known to play a key role in 
the genesis of abnormally rapid life-threatening heart rhythm 
disorders 1 1]. Despite considerable progress to date, complex 
spatiotemporal behaviors resulting from unstable spiral wave 
propagation remain poorlyunderstood theoretically with the 
exception of meander y, [3D, a classic spiral core instability 
with flower-like tip trajectories. A particularly rich dynam- 
ics results from instabilities in period-2 media where the local 
dynamics of the medium, i.e. the dynamics of uncoupled ex- 
citable elements, exhibits a period-doubling bifurcation as a 
function of parameters of the medium or the external stimula- 
tion frequency. Although period-2 behavior has been seen in 
different excitable and oscillatory media, it has received par- 
ticular attention in a cardiac context. The hallmark of period-2 
behavior in this context is alternans, a beat-to-beat alternation 
in the duration of cardiac excitation, which has been linked to 
the onset of lethal heart rhythm disorders |0]. 

Unstable spiral wave propagation in period-2 media is in- 
variably accompanied by "line-defects", which are the locus 
of points where the dynamics is locally period- 1 . Line-defects 
are generally present in these media when plane waves radi- 
ating out of the core region are unstable at the spiral rotation 
period, independently of whether meander is present or not. 
Studies in in vitro cardiac cell tissue cultures JJ, f^L chemical 
reactions 10, 0], and coupled oscillators [Oil 1M have re- 
vealed the existence of a rich variety of patterns ranging from 
one and three line-defect structures [0], to phase bubbles 
to line-defect turbulence lfl3ll . Spiral wave breakup in models 
of cardiac excitation has also been found in parameter regimes 
of local period-2 dynamics, and hypothesized in this context 
as a potential mechanism for heart fibrillation [0 [3, EH E3l • 
Spiral line-defect patterns, however, have not been systemati- 
cally investigated in cardiac tissue. 

In this Letter, we investigate the selection and dynamics of 
line-defect patterns resulting from unstable spiral wave prop- 
agation in cardiac tissue. Moreover, we interpret our findings 
using an amplitude equation framework recently used to study 
the evolution of line-defects during periodic stimulation from 
a single site fl7ll . In this framework, the spatiotemporal mod- 
ulation of the phase and amplitude of period-2 oscillations is 
described by a simple partial differential equation that can be 



readily analyzed. Our study is based on the standard wave 
equation for cardiac tissue 



d t V = jV 2 V-I m (V,y)/C n 



(1) 



where V is the transmembrane voltage, 7 is the voltage diffu- 
sion coefficient, C m is the membrane capacitance, and y is a 
vector of gate variables that controls the flow of ions through 
the membrane, and hence the total membrane ionic current 
I m . We studied different models of I m (V, y) and gating ki- 
netics to explore universal features of line-defect patterns that 
depend on qualitative properties of core and plane wave insta- 
bilities. The latter are manifested either as stationary 
or traveling W17I1 spatial modulations of period-2 oscillation 
amplitude with an intrinsic spatial scale determined by pa- 
rameters of the excitable medium W17I1 . These spatial modu- 
lations have nodes with period- 1 dynamics in one dimension, 
or nodal lines in two, which correspond here to line-defects 
in the spiral far-field. We therefore chose models to explore 
line-defect patterns for stationary and traveling nodes with and 
without meander. The model of Ref . fl5ll has pin wheel spirals 
(no meander) and stationary nodes under periodic pacing. The 
other two models of Ref. 11911 and Ref. 11711 both exhibit me- 
ander and have fixed and traveling nodes, respectively. 

Freely propagating spiral waves in all three models were 
studied by numerically solving Eq. (Q]) in a circular do- 
main of radius r e = 3 cm with no-flux boundary condi- 
tion, d r V\ r =r e = 0. Anchored spirals were studied by in- 
troducing an inexcitable disk of radius and imposing no- 
flux conditions on both the inner and outer radii, d r V\ r = ri — 
d r V\ r = r = 0. We implemented the phase-field method of 
Ref. B20ll that automatically handles no-flux boundary condi- 
tions in an arbitrary geometry using a finite-difference repre- 
sentation of the Laplacian on a square grid, and iterated Eq. (Q]) 
using a simple explicit Euler scheme. Model parameters are 
identical to the published ones except those listed in Fig. 1. 
The latter were chosen for intermediate action potential dura- 
tion restitution slopes, which suffice to produce unstable spi- 
ral waves with line-defects in each geometry, but are not steep 
enough to cause wave breakup in this domain size. 

We used a half plane wave as initial condition to initiate a 
spiral wave (obtained by first triggering a full plane wave and 
resetting part of the circular domain to the resting state). To 
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FIG. 1 : Membrane voltage (left) and corresponding alternans ampli- 
tude a (right) for different ionic models and geometries: model of 
Ref. fl5n without (a and b; Re = 1.15) and with a circular inex- 
citable disk of radius n = 0.5 cm (c and d; Re = 1.25), and model 
of Ref. fl9ll without disk (e and f; t/2 = 30 ms, t a % — 10 ms, 
th2 — 2.0 ms). In panels (a), (c), and (e), grey (orange online) indi- 
cates depolarization, and white indicates the trajectory of the spiral 
tip. In panels (b), (d), and (f), dark grey (red online) and light grey 
(yellow online) correspond to a > and a < regions, separated 
by line-defects (solid lines); the grey (color online) scale does not 
vary outside the indicated range of a. Line-defects in (b), (d), and 
(f), rotate counterclockwise with a period that is approximately 17, 
8, and 9 times the clockwise spiral rotation period in (a), (c), and 
(e), respectively. The discontinuity of a across the dashed line is a 
consequence of the definition of the common beat number (see text). 

track line defects, we define at each point x and time t a local 
beat number n(x, t), set everywhere initially to zero after the 
half plane wave is created, and increased by one at the end of 
each action potential, i.e. every time that the voltage V(x, t) 
crosses a fixed threshold V c with dV/dt < 0. We then define 
the period-2 alternans amplitude as 

a(x, t) = (-l) n ^ [D(x, n c (t)) - £>(x, n c (t) - 1)] /2 

(2) 

where Z>(x, n) = / v(x>t / )>Vc>n(X|t / )=n _ 1 df is the local ac- 
tion potential duration (APD) and n c (t) = min x n(x, t) is the 
common beat, i.e. the largest beat number that has been regis- 
tered at all points at time t. The line-defects are then the locus 
of points where a(x, t) = at any instant of time. The use 
of a common beat number introduces here a discontinuity in a 
(indicated by dashed lines in Fig. [T]) since the APD of a given 
beat might change as the wave front rotates around the spiral 




FIG. 2: Schematic illustration of transitions from (a) three to one 
and (b) one to three line-defects and corresponding examples from 
simulations of the model of Ref. (15] with the same conventions and 
parameters as in Fig.Q](b). The small circle (green online) in (a) and 
(b) represents the spiral core region or a small anchoring obstacle. 

tip. This discontinuity, however, does not affect the dynam- 
ics. Other methods to track line defects 10, [2]]] yield similar 
results except for unessential imaging differences. 

Results of simulations that pertain to the selection of the 
number of line-defects are shown in Fig. [I] The top four pan- 
els reveal that the pinwheel spirals simulated with the two- 
variable model of Ref. fl5ll exhibit three line-defects when 
propagating freely in spatially homogeneous tissue, but only 
one line-defect when anchored around an inexcitable disk of 
0.5 cm radius. Furthermore, the bottom two panels show that 
for the more physiologically realistic three-variable model of 
Ref. H, freely propagating spirals still exhibit three-line de- 
fects even though the spiral tip meanders. 

Since anchored spirals become free in the limit of vanish- 
ing obstacle size, one would expect transitions from one to 
three (three to one) line-defects to occur with decreasing (in- 
creasing) obstacle size. Indeed, for the model of Ref. W15I1 . we 
found three line defects for obstacles with diameter smaller 
than ~ 0.1 cm, including the freely propagating pinwheel spi- 
ral (ri = 0) in Fig. [T] (b), and one line defect for diameters 
larger than ~ 0.3 cm as in the example of Fig. [T] (d). For 
intermediate diameters, we found complex behaviors marked 
by transitions from three to one or one to three line-defects. 
The former occur when two line defects merge into one line 
defect that moves away from the core, and the latter when a 
phase bubble enclosed by a line-defect loop nucleates in, and 
expands from, the core, as illustrated in Fig. We find the 
same qualitative behavior in the model of Ref. L19] except 
that meander makes the transitions between patterns with dif- 
ferent numbers of line-defects more complex. 

Let us now turn to interpret our results in the amplitude 
equation framework 11711 . For simplicity, we restrict our anal- 
ysis to non-meandering spiral waves. Furthermore, to keep 
the analysis tractable, we first assume that the propagation 
wave speed is constant and relax this assumption subsequently 
when examining the motion of line defects. With this assump- 
tion, linear perturbations of a steady-state rigidly rotating spi- 
ral wave with period T obey the equation 

Td t a = aa + £ 2 V 2 a, (3) 

where a is the alternans amplitude subject to the radial 
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FIG. 3: Theoretical radial and angular (insets) dependence of the 
alternans amplitude for (a) 1 and (b) 3 line defect modes. 

d r a\ ri = da\ re = and angular a(6 + 27T, t) = —a(6,t) 
boundary conditions. The latter constrains the number of line 
defects to be an odd integer and results from the change in 
beat number across any closed circuit enclosing the spiral tip 
for steady- state alternans. It follows directly from the defini- 
tion of a [Eq. 0] and the requirement that the voltage be con- 
tinuous everywhere in space. In addition, a = In where 
the slope of the action potential duration restitution curve de- 
fined by D n+1 = f(T — D n ), controls the onset of alternans 
and £ ~ (7D) 1 / 2 , where D is the value of the action potential 
duration at the period-doubling bifurcation (a = 0), measures 
the scale over which the voltage dynamics is diffusively cou- 
pled on the time scale of one beat. 

This linear stability problem is easily solved by the sub- 
stitution a(r, t) ~ e m \I/(r, 6) that transforms Eq. © into a 
Helmholtz equation for \I/(r, 6). The latter can then be solved 
by separation of variables with the substitution \I>(r, 6) ~ 
R(r)Q(6). The angular part is found to be n (0) = 
sin ((n + 1/2)6), where mode n corresponds to 2n + 1 line- 
defects. The radial part obeys a Bessel equation. For vi > 0, it 
has solutions # n , m (r) oc J'_ n _ 1/2 (k 

n,m^ e) Jn-\-l/ 2 (^n,mX) 

J n+i/2( k n,rnr e )J-n-i/2(kn,mr) that satisfy the outer radial 
boundary condition da\ Te = 0, where n, m = 0,1,..., and 
the inner condition d r a\ ri =0 determines fc n>m , and hence the 
growth rate fl nym T = a — £,r~ 2 k 2 m . We find that the small- 
est fc n>m occurs for n = independently of the ratio r e /ri. 
Therefore, the mode corresponding to a single line defect is 
the most unstable when the spiral is anchored. This agrees 
with our numerical observations in Fig.[T](d). For freely rotat- 
ing spirals, T{ = 0, J_ n _i/ 2 (r) diverges at the origin, so the 
solutions are R n ,m(r) oc J n+1 / 2 (k nj7n r), where k nm r e is the 
m th zero of J' n (r). The most unstable modes are n = and 
n = 1 corresponding to 1 and 3 line defects, respectively (see 
Fig. [3]). However, Ji 72(^0,0^) has a divergent derivative that 
is incompatible with the physical requirement that the volt- 
age, and hence the APD, must vary smoothly on a scale £. On 
the other hand, ^3/2(^1,0^) smoothly vanishes at the origin. 
Therefore, in this case, the boundary condition at the origin 
selects a 3 -line-defect pattern as observed in Fig.[T](b). Inter- 
estingly, a 3 -line-defect pattern is also selected with meander 
present [Fig.[T](f)], thereby suggesting that the boundary con- 
dition on a on the outer scale of the line-defect pattern is not 
strongly affected by meander. 

The analysis also predicts qualitative features of the radial 
distribution of alternans amplitude for three- and one-line- 
defect patterns of Figs. [TJb) and Old), respectively. Fig. [4] 
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FIG. 4: (a)rms versus radial distance r for the pinwheel spiral (thin 
solid line, red online) and anchored spiral (thick solid line, blue on- 
line), and respective theoretical radial modes Ri,o(r) (thin dashed, 
red online), and Ro,o(r) (dotted-dashed line, blue online), 
compares the numerical radial distributions of root-mean- 
square amplitude (a) rms averaged over a full line-defect rota- 
tion period for three line-defects (thin solid line) and one line- 
defect (thick line), with the corresponding radial modes from 
the theoretical analysis (dashed lines) scaled to have the same 
radial average as the observed curves. The theory predicts 
well that the alternans amplitude is more strongly suppressed 
near the core for the larger number of line defects. 

So far our analysis has assumed that the wave speed is con- 
stant, which predicts that line-defects extend straight out of 
the core and are stationary, as implied by the angular distribu- 
tion sin ((n + 1/2)6) of linearly unstable modes (see Fig. [3]). 
In contrast, simulations in Fig. [I] show that line-defects have a 
spiral shape and slowly rotate inward in the opposite direction 
of the spiral wavefront. Line-defect motion can generally be 
induced both by line-defect curvature and the dependence of 
the wave speed c on the interval / between two waves, known 
as the conduction velocity (CV) restitution curve in the car- 
diac literature. While a full stability analysis that includes 
these effects would be required to treat line-defect motion in 
general, two important limiting cases can be readily analyzed. 

The first pertains to anchored spiral waves for medium pa- 
rameters where plane waves paced at the spiral rotation period 
exhibits stationary line-defects, as for the model of Ref. fl5ll 
studied here. In this case, we expect line-defect motion to 
be generated predominantly by the spiral wavefront dynamics 
around the anchoring obstacle. Neglecting wavefront curva- 
ture effects, this dynamics should be approximately described 
by that of a propagating pulse in a one-dimensional ring of 
perimeter L = 2nri uj\, l22|]. To test this hypothesis, we 
computed the quasiperiodic frequency Q of the local medium 
dynamics induced by line-defect rotation for anchored spirals 
for the model of Ref. fl5ll . The frequency was obtained 
by fitting the time series a(r, jT)/a(r, 0) at a single point 
r to rf cos(QTj + S), with 77, Q, and S the fitting parame- 
ters. For the theory, we used the dispersion relation giving the 
quasiperiodic frequency ft modulating alternans, a ex e*^ jT , 
in a one dimensional ring derived in Ref. iTTtH 
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FIG. 5: Comparison of numerical (circles, red online) and theoretical 
(solid line) line-defect rotation frequencies for spiral waves anchored 
to disks of varying radii for the model of Ref . II 1511 . 
where k = tt/L + flT/L is the wavenumber corresponding 
to a single line defect and A = c / (/)/(2c 2 ). The APD- and 
CV-restitution curves, /(/) and c(I), were calculated in a one- 
dimensional cable as in Ref . 117]. In addition, the intercellular 
coupling parameters uo and £ were estimated as uo ~ Zj/c 
and£ ~ (7^) 1/2 & 

The comparison in Fig. \5\ shows that 
the ring-based theory predicts reasonably well the frequency 
of line-defect rotation for anchored spiral waves of different 
period T, which was varied here by increasing the obstacle 
radius vi in the simulations. 

The opposite limit that can also be readily understood is 
the one where plane waves paced at the spiral rotation period 
exhibit line-defects that move towards the pacing site, which 
generally occurs for steeper CV-restitution. In this case, line- 
defect motion is expected to be dominated by the far-field spi- 
ral dynamics 111711. We have checked that, for the two-variable 
model of Ref. [17], spiral line-defects indeed rotate inward 
with a frequency equal to the product of the velocity of the 
planar line-defects and the inverse of their spacing. This prop- 
erty was purposely checked in a domain much larger than the 
spiral wavelength (r e = 18 cm) and with an obstacle size 
(Xi = 0.72 cm) sufficient to prevent spiral wave breakup in- 
herent in this model. However, we expect this behavior to be 
generic for systems with traveling planar line-defects and to 
also apply to freely rotating spirals with three line defects for 
parameters where breakup does not occur. 

In summary, we have surveyed spiral line-defect patterns 
in simplified models of cardiac excitation with period-2 dy- 
namics. Although far from exhaustive, this survey yields the 
striking finding that freely propagating and anchored spiral 
waves select different numbers of line-defects. This opens up 
the possibility to distinguish free and anchored spiral waves 
in cardiac tissue by monitoring the number of line-defects. 
We have shown that spiral wave unstable modes with differ- 
ent numbers of line-defects correspond to topologically quan- 
tized solutions of a Helmholtz equation. In this framework, 
the boundary condition on the period-2 oscillation amplitude 
in the spiral core, which is fundamentally different for free 



and anchored spirals, is responsible for selecting the num- 
ber of line defects. Furthermore, we have found that spiral 
line-defect inward rotation can be driven either by the core 
or far- field wavefront dynamics, with concomitantly differ- 
ent frequencies. Our results suggest that the observation of 
single-line-defect spirals in cardiac tissue culture may 
be a consequence of anchoring on small millimeter- size het- 
erogeneities. However, the dynamics in real tissue is also in- 
fluenced by the coupling of voltage and intracellular calcium 
dynamics lll|4j], which has been neglected here. The investi- 
gation of the effect of this coupling on line-defect dynamics 
and its relationship to wave breakup is an interesting future 
project. Finally, the previous finding of free spirals with one 
line-defect 112 ill in period-2 media with qualitatively different 
excitable dynamics than cardiac tissue suggests that other pat- 
tern selection mechanisms may be operative in different me- 
dia. These differences also remain to be elucidated. 

We thank Bias Echebarria for valuable discussions. This 
work was supported by NIH Grant No. P01 HL078931. 
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